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Abstract 



A search for signatures of extra dimensions in the diphoton invariant-mass spectrum 
has been performed with the CMS detector at the LHC. No excess of events above 
the standard model expectation is observed using a data sample collected in proton- 
proton collisions at a/s = 7TeV corresponding to an integrated luminosity of 2.2 fb^^. 
In the context of the large-extra-dimensions model, lower limits are set on the effective 
Planck scale in the range of 2.3-3.8 TeV at the 95% confidence level. These limits are 
the most restrictive bounds on virtual-graviton exchange to date. The most restrictive 
lower limits to date are also set on the mass of the first graviton excitation in the 
Randall-Sundrum model in the range of 0.86-1.84 TeV, for values of the associated 
coupling parameter between 0.01 and 0.10. 
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Over a decade ago, Arkani-Hamed, Dimopoulos, and Dvali (ADD) HI |2l proposed that ex- 
tra spatial dimensions could potentially solve the hierarchy problem in the standard model 
(SM) [3J. They proposed a scenario in which the SM is constrained to the common 3+1 space- 
time dimensions (brane), while gravity is free to propagate throughout a larger multidimen- 
sional space (bulk). The gravitational flux on the brane is therefore diluted by virtue of Gauss's 
Law in the bulk, which relates the fundamental Planck scale Md to the apparent reduced scale 
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Mpi ^2x 10^^ GeV according to the formula Md"™^^ = where r and Med are the size and 
number of the extra dimensions (ED), respectively. Postulating a fundamental Planck scale 
to be on the order of the electro weak symmetry breaking scale (~1 TeV) results in ED with 
r < 1 mm for Hed > 2. 

Another model of ED that solves the hierarchy problem was proposed by Randall and Sun- 
drum (RS) m. In this model, the observed hierarchy is related instead to the warped anti-de 
Sitter (AdS) geometry of the ED. We specifically consider the RSI model whereby only one fi- 
nite ED exists separating two branes, one at each end of the ED. The geometry of the bulk is 
based on a slice of a five-dimensional AdS space with a length ttTc, where is the compactifi- 
cation radius. The full metric is given by ds^ = e^^^^Vfj^y dx^dx^ — r'^dy^, where Greek indices 
run over 4-dimensional space-time, r]^y is the Minkowski metric tensor, and < i/ < tt is the 
coordinate along the single ED of radius r^-. The value of k specifies the curvature scale (or 
"warp factor ") and relates the fundamental Planck scale on one brane to the apparent scale on 
the other by A^^ = M^ie"^^"'^ . As a consequence, a value of krc ~ 10 would provide a natural 
solution to the hierarchy problem, yielding A^^ ~ 1 TeV. 

Phenomenologically, the excited gravitons in both models preferentially decay into two gauge 
bosons, such as photons, rather than into two leptons, because the graviton has spin 2, and 
so fermions cannot be produced in an s wave. In the RS scenario, gravitons appear as well- 
separated Kaluza-Klein (KK) excitations with masses and widths determined by the parame- 
ters of the RSI model. One convenient choice of parameterization is the mass Mi of the first 
excitation of the graviton and the dimensionless warp factor k = k/M^i, which defines the 
strength of associated coupling to the SM fields. Precision electroweak data constrain k > 0.01, 
while perturbativity requirements limit k < 0.10. 

In the ADD model, the wavefunction of the KK gravitons must satisfy periodic boimdary con- 
ditions, resulting in discrete energy levels with modal spacing of the order of the inverse ED 
size, from 1 meV to 100 MeV, much smaller than the spacing in the RSI model, which is ex- 
pected to be of the order of 1 TeV. This effect produces an apparent continuum spectrum of 
diphotons, rather than distinct resonances, at high (~1 TeV) diphoton invariant mass M^^. 

Summing over all KK modes in the ADD scenario results in a divergence in the cross section, so 
an ultraviolet (UV) cutoff scale Ms is imposed. This effective Planck scale is related to — ^but po- 
tentially different from — the fundamental Planck scale Md- The precise relationship depends 
on the UV completion of the effective theory. The effects of virtual-graviton production on the 
differential diphoton cross section are parameterized by the single variable rjQ = J-"/ Mg, where 
jF is an order-unity dimensionless parameter, for which several conventions exist 15] . In this 
Letter, we measure lower limits on Ms in three different conventions: GRW [6|, Hewett fT\, and 
HLZ |8|. 

Searches for ED via virtual-graviton effects in the ADD model have been conducted at HERA, 
LEP, the Tevatron, and the LHC f9j[l0l. The most stringent previously published limits on Ms 
for Wed > 3 come from the previous measurement in the diphoton channel at the Compact 
Muon Solenoid (CMS) experiment [5J. For 

f^ED — 2, the DO measurements in the dijet |fTl| and 
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diphoton+dielectron ||T2| channels are more restrictive. The most sensitive previous search for 
RS gravitons was conducted by the ATLAS experiment ||l3l . They used a search in the dilepton 
final state to exclude Mi < 1.63 TeV for k = 0.10. 

In this Letter, we present a search for both non-resonant and resonant diphoton production, in 
the ADD and RSI models, respectively. We use data corresponding to an integrated luminosity 
of 2.2 fb^^, collected in pp collisions at = 7 TeV at the LHC with the CMS detector between 
March and August 2011. 

The CMS detector llMll is designed to study collisions at the LHC. An all-silicon tracker, an 
electromagnetic calorimeter (ECAL), and a hadronic sampling calorimeter (HCAL) are all con- 
tained within a large-bore 3.8 T superconducting solenoid. In the central region, the tracker 
consists of three radial layers of silicon pixel detectors followed radially by silicon strip detec- 
tors. The finely segmented ECAL has a design resolution for unconverted photons better than 
0.5% at energies exceeding 100 GeV in the barrel {\rj\ < 1.44). Here, the pseudorapidity rj is 
defined as — ln(tan(0/2)), where 6 is the polar angle with respect to the direction of the coun- 
terclockwise beam. Beyond the solenoid lie four layers of muon detectors. The instantaneous 
luminosity is measured with a relative uncertainty of 4.5% using information from forward 
hadronic calorimeters fT5]. The two-tiered CMS trigger consists of the level-one trigger, com- 
posed of custom hardware, and the software-based high-level trigger. 

Events for the analysis were collected through a diphoton trigger, where each photon was re- 
quired to have a transverse energy Ej = E sin of at least 33, 55, or 60 GeV, depending on 
the instantaneous luminosity. We require that an event be consistent with a pp collision and 
have at least one well-reconstructed primary vertex |1T6| . We then reconstruct photons with 
Ej > 70 GeV in the ECAL barrel by clustering electromagnetic energy depositions. Electrons 
that do not originate from photon conversions are suppressed by using information from the 
pixel detector to associate tracks and ECAL clusters compatible with an electron hypothesis. 

Hadronic jets can be misidentified as photons when their leading hadron is an energetic or 
rj meson. We reduce the misidentification rate from this source by placing the same restric- 
tions on the isolation as in the previous analysis for this channel fS]. These restrictions limit 
the total transverse energy because of tracks and calorimeter depositions near the photon clus- 
ter. Restrictions on the shower-shape variable CTtjrj, which is a modified second moment of the 
electromagnetic energy cluster about its mean rj position [T7J, also suppress hadronic misiden- 
tification. Topological and timing criteria suppress anomalous signals present in the ECAL [18J. 
Diphoton events are selected in which My-y > 140 GeV. 

The photon reconstruction and identification efficiency is determined in Monte Carlo (MC) sim- 
ulation and corrected using a data-to-MC scale factor derived from studying Z — > e+e^ events. 
The measured efficiency for a single Ej > 70 GeV photon with \t]\ < 1.44 is (87.4 ± 5.4)% and 
depends only weakly on the Ej and t] of the photon, and the number of extra collisions present 
in the event. The systematic uncertainty covers the variation as a function of these variables, 
the most significant of which is the number of extra collisions. We reweight the simulation to 
give the same reconstructed primary-vertex distribution (on average 6-8 vertices) as observed 
in the data. We determine the corresponding diphoton reconstruction and identification effi- 
ciency (76.4 ± 9.6)% by squaring the single-photon efficiency. 

The simulation of ED in the ADD model is performed using version 1.3.0 of the SHERPA ||T9| 
MC generator. The simulation includes both SM diphoton production and signal diphoton pro- 
duction via virtual-graviton exchange in order to account for the interference effects between 
the SM and ADD processes. The leading-order (LO) SHERPA cross sections are multiplied by a 
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constant next-to-leading-order (NLO) K factor of 1.6 ± 0.1, a value that represents an updated 
calculation by the authors of Refs. [1201 121 1. The systematic uncertainty on the signal K factor 
reflects the approximate variation of the K factor over a large region of the model parameters; 
it is not intended to account for the theoretical uncertainty. The cross sections in the simulation 
are conservatively set to zero for > Ms because the theory becomes non-perturbative for 
larger values of V s. Introducing this sharp truncation reduces the upper limits on Ms by a few 
percent. 

The simulation of RS-graviton production is performed using version 6.424 of the PYTHIA ||22| 
MC program. The signal cross section is scaled by a mass-dependent NLO K factor Il20l 121} , 
which ranges from 1.6 to 1.8 as a function of Myy and for different values of k. The CTEQ6l1 Il23l 
parton distribution functions (PDF) are used in the simulation of both the ADD and RS models, 
and a 1.5% relative uncertainty on the signal acceptance is included, associated with the choice 
of PDF. 

Optimization of the event selection is done separately for both ADD and RS scenarios. The 
signal in both cases is predominantly at central values of rj, while the high-M-y-y SM background 
dominates the signal in the forward region; therefore, we restrict ourselves to photons located 
in the ECAL barrel only. In the ADD scenario, we find that the optimal requirement on M-y-y, 
based on the discovery potential, is 900 GeV. 

In the search for RS gravitons, a fixed window is selected about the Mi mass point of interest. 
Because the signal shapes deviate from Gaussian distributions, we define an effective measure 
of the signal width c^gff as the half-width of the narrowest mass interval containing 68% of the 
signal from simulation. A window is then formed about the resonance mean of size ±5 a^ff in 
the data. This window contains 96-97% of the signal acceptance for all mass points considered 
in this analysis, and the detector resolution is negligible with respect to the window size. This 
choice of the window maximizes the signal acceptance and analysis sensitivity in the case of 
small backgrounds. 

Backgrounds from the misidentification of a hadronic jet as a photon are small in the signal 
region but contribute to the low-M-y-y region. Two such sources of backgrounds from isolated- 
photon misidentification are considered: multijet production and prompt single-photon (7 + 
jet) production. In particular, we measure on a background-dominated sample a misidentifi- 
cation rate, defined as the ratio of the number of isolated photon candidates to non-isolated 
photon-like objects. These photon-like objects are reconstructed as photons but fail one of the 
isolation or shower-shape criteria; therefore, the samples corresponding to numerator and de- 
nominator are mutually exclusive, and prompt photons have a negligibly small contribution 
to the denominator. The misidentification rate is measured in a photon-triggered sample in 
bins of photon(-like) candidate Ej, but the objects used in the measurement are required to be 
well-separated from the triggered object to avoid a trigger-induced bias. A cross-check of the 
misidentification rate is performed in a jet-triggered sample. 

Because the background-dominated sample in which we measure the misidentification rate 
may contain some genuine, isolated photons that "contaminate" the numerator of the misiden- 
tification rate, we correct for this on a bin-by-bin basis. The cr^,j requirement is released and the 
numerator sample is fit for the fraction of prompt photons using one-dimensional probability 
density histograms ("templates") in CTjjij. The signal template is constructed from MC simula- 
tion, and the background template is constructed from reconstructed photons that fail one or 
more of the isolation criteria. The measured misidentification rate falls from 7% at Ej = 70 
GeV to 2% at Ej = 120 GeV. We apply a 20% systematic uncertainty to the rate derived from 
the variation of the misidentification rate measured in a jet-triggered sample. 
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The multijet and 7 + jet backgrounds to the reconstructed diphoton spectrum are estimated 
by using the misidentification rate to extrapolate from two background-dominated reference 
regions, both selected with the same diphoton trigger as the primary signal sample. One re- 
gion includes events with only one isolated photon, but one or more non-isolated photons. The 
other region includes events with no isolated photons, but two or more non-isolated photons. 
The diphoton trigger is sufficiently inclusive that the regions are unaffected by the trigger se- 
lection. By applying the prompt-photon misidentification rate to these two reference regions, 
we predict the 7 + jet and multijet backgrounds in the signal region. 

The SM diphoton background dominates the signal region. The expected number of back- 
ground events due to this process is computed by rescaling the prediction from pythia with 
an NLO K factor that varies with M^-y. The NLO prediction is calculated with the DIPHOX+ 
GAMMA2mc Il24ll25| generators, which take into account the fragmentation processes in which 
the photons can come from the collinear fragmentations of hard partons. A separate analy- 
sis by CMS has also demonstrated good agreement with the NLO prediction at low M-y-y < 
300 GeV [26|. The sub-leading-order gluon-fusion box diagram is included as a part of the 
PYTHIA calculation because of its large contribution at the LHC energy, although its effects are 
small at high M-y-y. The K factor varies between 1.7 and 1.1 from low to high M-y-y. A sys- 
tematic uncertainty of 15% on the value of the K factor is determined by examining the PDF 
uncertainties and variation of the renormalization and factorization scales. 




200 400 600 800 1000 1200 1400 1600 1800 2000 
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Figure 1: Observed event yields (points with error bars) and background expectations (filled 
solid histograms) as a function of the diphoton invariant mass. Photons are required to be 
isolated, with £t > 70 GeV and < 1.44. The shaded band around the background estimation 
corresponds to the systematic uncertainty. The last bin includes the sum of all contributions 
for Myy > 2.0 TeV. The simulated distributions for two, non-excluded signal hypotheses are 
shown for comparison as dotted and dashed lines. 

Figure [l] shows the invariant mass distribution of the selected events, together with the es- 
timated distributions for each of the backgrounds. Table [l] presents the observed number of 
events in the data and the predicted number of background events in different ranges in Myy 
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and corresponds directly to Fig. [T] The last column corresponds to the signal region for the 
ADD search. We find that the observed data are consistent with the background estimate 
throughout the M^^ spectrum and do not show an excess of events, neither resonant nor non- 
resonant. 

Table 1: Observed event yields and background expectations for different reconstructed dipho- 
ton invariant-mass ranges. Full systematic uncertainties are included. 



Process 




Diphoton Invariant Mass Range [TeV] 




[0.14, 0.2] 


[0.2, 0.5] 


m 


h 0.9] 




0. 


9 


Multijet 


15 


±6 


17 ± 7 


0.2 


±0.1 


0.003 


± 


0.001 


7 + jet 


102 


± 15 


124 ± 18 


2.5 


±0.4 


0.19 


± 


0.04 


Diphoton 


372 


± 70 


414 ± 78 


16.9 


± 3.2 


1.3 


± 


0.3 


Backgrounds 


489 


± 73 


555 ±81 


19.6 


± 3.2 


1.5 


± 


0.3 


Observed 


484 


517 




16 




2 





To set limits on virtual-graviton exchange in the ADD scenario, we compare the number of 
observed and expected events in the signal region (M-y-y > 0.9 TeV) and set 95% confidence 
level (CL) upper limits on the quantity S = (c^totai ~ <^sm) x B x A, where dtotai represents the 
total diphoton production cross section (including signal, SM, and interference effects), and 
(TgM represents the SM diphoton production cross section. The signal branching fraction to 
diphotons is indicated by B and the signal acceptance by A. We use the CLg technique to 
compute the limits with a likelihood constructed from the Poisson probability to observe N 
events, given S, the signal efficiency (76.4 ± 9.6)%, the expected number of background events 
(1.5 ± 0.3), and the integrated luminosity C = (2.2 ± 0.1) fb"^ lilSj. 
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Figure 2: Signal cross section S parameterization as a function of the strength of the ED effects, 
rjG (left) and as a function of 1/Mg for the HLZ = 2 case (right). 

The observed 95% CL upper limit on S is 3.0 fb. For the HLZ Hed = 2 case, we parameterize 
S directly as a smooth function of 1/Mg. For all other conventions, S is parameterized as a 
function of the parameter rjc, as in Ref. |5|. The observed 95% CL limit, together with the 
signal parameterization, are shown in Fig. |2] The intersection of the cross-section limit with 
the parameterized curve determines the 95% CL upper limit on the parameter r]Q. As seen 
from the plot, these upper limits on S correspond to upper limits of f/c < 0.0097 TeV^^ and 
1/M| < 0.0055 TeV"^. The upper limits on rjc are equated to lower limits on Mg and are 
shown in Table |2l 
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Table 2: The 95% CL lower limits on Ms (in TeV) in the GRW, Hewett, and HLZ conventions for 
two values of the ADD signal K factor. All limits are computed with a signal cross section trun- 
cated to zero for \/s > Mg. The limits are presented for both positive and negative interference 
in the Hewett convention and for rivu = 2-7 in the HLZ convention. 



K factor 


GRW 


Hewett 
pos. neg. 


2 


3 


HLZ (med) 
4 5 


6 


7 


1.0 


2.94 


2.63 2.28 


3.29 


3.50 


2.94 2.66 


2.47 


2.34 


1.6 


3.18 


2.84 2.41 


3.68 


3.79 


3.18 2.88 


2.68 


2.53 




[GeV] 

Figure 3: The 95% CL exclusion region for the RSI graviton model in the M^-k plane. The 
expected limits coincide very closely with the measured limits and so are not shown in the 
figure. Also shown are bounds due to electroweak constraints l,27J and naturalness {An > 
10 TeV). Perturbativity requirements bound k < 0.10. 
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For the RS scenario, the same limit-setting calculation is performed, but in a bounded window 
in Myy. The expected upper limits coincide within a few GeV of the measured upper limits. 
Figure [3] shows the excluded regions in the Mi-k plane. Also shown are bounds due to preci- 
sion electroweak measurements [27J and to naturalness arguments. Table |3] presents the 95% 
CL lower limits on the graviton mass Mi for different values of k. 

Table 3: The 95% CL lower limits on Mi for given values of the coupling parameter, k. For 
k < 0.03, masses above the presented limits are excluded by electroweak and naturalness con- 
straints. 



k 


0.01 


0.02 


0.03 


0.04 


0.05 


0.06 


0.07 


0.08 


0.09 


0.10 0.11 


Ml [TeV] 


0.86 


1.13 


1.27 


1.39 


1.50 


1.59 


1.67 


1.74 


1.80 


1.84 1.88 



In summary, we have performed a search for extra spatial dimensions leading to enhanced 
resonant or non-resonant diphoton production in proton-proton collisions at a center-of-mass 
energy of 7 TeV at the LHC. Using a data sample corresponding to an integrated luminosity of 
2.2 fb^^ recorded by the CMS experiment, we observe no excess in diphoton production above 
the rate predicted from SM background sources. Values of the effective Planck scale Ms less 
than 2.3-3.8 TeV are excluded at 95% CL for ADD models. We also exclude at 95% CL resonant 
graviton production in the RSI model with values of Mi less than 0.86-1.84 TeV depending on 
the normalized coupling strength, k. We present limits on both the ADD and RSI models of 
extra dimensions in the diphoton final state that extend those observed at the Tevatron |[T2|, as 
well as those set previously by the CMS [5 1 and ATLAS IIT3II experiments. 
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